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p/j ' Abstract 

In this paper we generalize and analyze the model for pricing American-style 
Asian options due to (Hansen and J0rgensen 2000) by including a continuous divi- 
i" 1 dend rate q and a general method of averaging of the floating strike. We focus on the 

qualitative and quantitative analysis of the early exercise boundary. The first order 
Taylor series expansion of the early exercise boundary close to expiry is constructed. 
We furthermore propose an efficient numerical algorithm for determining the early 
exercise boundary position based on the front fixing method. Construction of the 
£S| | algorithm is based on a solution to a nonlocal parabolic partial differential equation 

for the transformed variable representing the synthesized portfolio. Various numer- 
ical results and comparisons of our numerical method and the method developed by 
(Dai and Kwok 2006) are presented. 
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1 Introduction 

Evolution of trading systems influences the development of financial derivatives market. 
First, simple derivatives like forwards and plain vanilla options were used to hedge the 
risk of a portfolio. Progress in pricing these simple financial instruments pushed traders 
into inventing less predictable and more complex derivatives. Using financial derivatives 
with more complicated pay-offs brings into attention also new mathematical problems. 
Asian options belong to a group of the so-called path-dependent options. Their pay-off 
diagrams depend on the spot value of the underlying during the whole or some part(s) of 
the life span of the option. Usually Asian options depend on the (arithmetic or geometric) 
average of the spot price of the underlying. They can be also used as a useful tool for 
hedging highly volatile assets or goods. Since the price of an underlying varies during the 
life span of the option the holder of the Asian option can be secured for the risk a sudden 
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price jumping to an undesirable region (too high for the call option holder or too low for 
the put option holder). Among path dependent options Asian options plays an important 
role as they are quite common in currency and commodity markets like e.g. oil industry 
(cf. (Wilmott, Howison and Dewynne 1995, Hull 1997, Wu, Kwok and Yu 1999, Hansen 
and J0rgensen 2000, Detemple 2006, Dai and Kwok 2006, Wystup 2006, Kwok 2008, Kim 
and Oh 2004, Wu and Fu 2003, Linetsky 2004)). 

In this paper we focus on the so-called floating strike Asian call or put options whose 
strike price depends on the averaged path history of the underlying asset. More precisely, 
we are interested in pricing American-style Asian call and put options having the pay-off 
functions Vt(S, A) = (S — A) + and Vr(S, A) = (A — S) + , resp. The strike price A is given 
as an average of the underlying over the time history [0, T\. We are analyzing the early 
exercise boundary for American-style Asian option (cf. (Hansen and J0rgensen 2000, Dai 
and Kwok 2006, Kwok 2008)). Recall that American-style options can be exercised at 
any time until the maturity T. The holder of such an option has the right to exercise it 
or to keep it depending on the spot price of the underlying S t at time t and its history 
{5'm ; < u < t} prior the time t. The boundary between "continuation" and "stopping" 
regions plays an important role in pricing American-style of options. It can be described 
by the mapping : t i— > S**, where S% is the so-called early exercise boundary (cf. 
(Hull 1997, Geske and Johnson 1984, Geske and Roll 1984, Karatzas 1988, Chadam 2008, 
Kwok 2008, Kuske and Keller 1998, Mallier 2002, Pascucci 2008). 

The paper is organized as follows. In the next section, we discuss the probabilistic 
model for valuation of American type Asian options with a floating strike given in the form 
of an average of the underlying asset price. The model is based on conditioned expected 
values and theory of martingales. In the third section, we present a general result enabling 
us to calculate the value of the limit of the early exercise boundary at expiry. We also 
calculate the analytical integral formula for an option with continuous geometric average 
and an approximation for the value of an option with a continuous arithmetic average. 
The main result of this section is the approximation formula for the first order Taylor 
series approximation of the early exercise boundary near to expiry. Similarly as in the 
case of plain vanilla call options, we show that the leading order of the expansion is the 
square root \JT — t of the time T — t remaining to maturity. In the fourth section, we 
present the fixed domain transformation method yielding a nonlocal nonlinear parabolic 
partial differential equation for pricing the synthesized portfolio for an Asian option. We 
also present an efficient and robust numerical scheme for construction of an approximation 
of the solution to the governing system consisting of a partial differential equation with 
an algebraic constraint. In the last part of this section we present also the results of the 
presented method on some examples. 

2 A probabilistic approach for pricing of American- 
Style of Asian options 

The main purpose of this section is to derive an integral equation for valuation of the 
early exercise boundary of an American-style Asian option paying continuous dividends. 
We follow the ideas of derivation due to (Hansen and J0rgensen 2000). Their formula 
for a floating strike option was derived using the theory of martingales and conditioned 
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expected values. We extend their formula to Asian options on underlying paying non-zero 
dividend yield and having a general form of floating strike averaging. In a more detail, 
we discuss geometric, arithmetic and weighted arithmetic averaging operator. 

The pricing model is based on the assumption on the stochastic behavior of the un- 
derlying asset in time. Throughout the paper we shall assume that the underlying asset 
price St is driven by a stochastic process satisfying the following stochastic differential 
equation 

dS t = (r - q)S t dt + aS t dW t p , < t < T. (1) 

It starts almost surely from the initial price Sq > 0. Here the constant parameter r > 
denotes the risk-free interest rate whereas q > is a continuous dividend rate. The 
constant parameter a stands for the volatility of the underlying asset returns and Wf is a 
standard Brownian motion with respect to the standard risk-neutral probability measure 
P. A solution to equation ([1]) corresponds to the geometric Brownian motion 

S t = S e( r - q -^ t+aW t P , 0<t<T. 

We shall derive an integral equation determining the value of an American-style Asian 
option with a floating strike. If we define the optimal stopping time as T*, the pay-off of 
the option is set by 

V T *= (p(5 T . -A t *)) + , 

where Vt is the value of the option at time t, A t is a continuous average of the underlying 
asset value during the interval [0, t] and p = 1 for a call option and p = —1 for a put option. 
We may consider several different types of continuous averages presented in Table [TJ 



Table 1: Classification of averaging methods. 



arithmetic average 


geometric average 


weighted arithmetic average 


A t = \ f* S u du 


In At = ~ /q In S u du 


At ~ fiats) dsfo^ U ) S - du 



In the case of a weighted arithmetically averaged floating strike Asian option the kernel 
function a(.) > with the property J °° a(s) ds < oo is usually defined as a(s) = e~ Xs 
where A > is constant. 

According to (Hansen and J0rgensen 2000), the American-style contingent claims can 
be priced by the conditioned expectations approach. The option price can be calculated 
by considering all possible stopping times in the interval [t, T] 



S t = S,A t = A 



V(t, S, A) = esssup Ef [e"^^ (p(S s - A s 

where 7[t,T\ denotes the set of all stopping times in the interval [t,T] and Ef[X] = 
E p [A|jF 4 ] is the conditioned expectation with information available at time t (the in- 
formation set is represented by the filtration Tt of the cr-algebra T where the Brownian 
motion is supported). To simplify the formula we change the probability measure by the 
martingale 

So 
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the new probability measure Q being defined as dQ = rixdP. According to Girsanov's 
theorem (Revuz and Yor 2005), the process 



W t p - at 



is a standard Brownian motion with respect to the measure Q. The value of the underlying 
asset price under this measure is defined by 



S t = 5ne (^+^ 2 )*+^ s . 



(2) 



All assets priced under this measure are Q-martingales when discounted by the underlying 
price. According to this fact, we can reduce the dimension of stochastic variables. We 
introduce a new variable xt 



j±. We obtain: 



V(t,S,A) 



essup Ef e - r(s -' } [p(S s - A s 
esssup E® \2L e -r{»-t) (p(S 8 - A, 



It.T] 



esssup Ej 

SGT lt,T] 



Q. 



esssup e^^SEf (p(l 



Xq 



S t = S,A t = A 

S t = S,A t = A 
S t = S,A t = A 
S t = S,A t = A 



The last expression can be rewritten in terms of the new variable X ~q CIS follows: 



V(t,x) 



-<jt 



V(t, S, A) _ T * 



S 



E; 



p(l - x T * 



where T 4 * = inf{s e [t,T]\x s = x*} and the function [0,T] 3 t h-> x* t G 
early exercise boundary. 



(3) 

describes the 



Definition 2.1 The stopping region S and continuation region C for American- style 
Asian call and put options (T3|) are defined by 

Scan = C put = {(t,x)\t E [0,T],0 < x < x* t }, Caoi = S put = {{t,x)\t E [0,T),x* < x < oo}. 

where [0, T] G 1 1— > x* t G R zs a continuous junction determining the early exercise bound- 
ary. By ls(-) we s/ia/Z denote the indicator function of the set S, i.e. ls{t,y) = 1 for 
(t,y) G S and ls(t,y) = otherwise. 

In the following theorem we present a solution to the pricing problem with one stochas- 
tic variable Xt formulated in Q. It is a generalization of the result by (Hansen and 
J0rgensen 2000) and (Wu et al. 1999) for the case of a nontrivial dividend rate q > and 
a general form of the averaging of the floating strike price. 
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Theorem 2.1 The value V(t, x t ) of the American-style floating strike Asian call (p = 1) 
or put option (p = — 1) on the underlying asset Xt paying continuous dividends with a rate 
q > is given by V(t, x t ) = v(t, x t ) + e(t, x t ), where 



v(t,x t ) 
e(t, x t ) 



E; 



E; 



pe~ qu x„ 



u l s (u, x u ) - (r - qx u l )duj 



(4) 
(5) 



with the average given by the function A t and stopping region S. 
In the proof of Theorem 12.11 we shall use the next lemma. 

Lemma 2.1 The auxiliary variable Xt = 4r satisfies the following stochastic differential 
equation: 



dx t 



dA f 



Xf 



At 



[r — q)x t dt — ax t dW^. 



Proof: [of Lemma 12. lj We express the differential dx t = d 
1 . . A t ,„ , A t . m , 2 dA t 



as 



dxt 



5, 



■dA t - '-^dSt + '-^(dSt) 2 



Xf 



— (r — q)xt dt — ox t dW®, 



and the proof of lemma follows. □ 

Notice that, when comparing to the original expression due to (Hansen and J0rgensen 
2000) with a zero dividend rate q = the only difference is that the parameter r is 
replaced by the term r — q. 

The value of ^ depends on the method of averaging of the underlying asset used 
in valuation. The expressions for continuous averages are presented in Table |2j In this 
table we present the value of the weighted arithmetic average restricted to the kernel 



as 



-As 



Table 2: The value of the differential ^r*. 

At 



arithmetic average 


geometric average 


weighted arithmetic average 




dA t 1 1 9 j-i 

—rn- = — -t In Xt at 

A t * * 





Proof: [of Theorem 12. lj We follow the proof of the original result by (Hansen and 
J0rgensen 2000) and we include necessary modifications related to the form of averaging 
and the fact that q > 0. 

First, we suppose that (t, x) G C. The option is held and so we can apply Ito's lemma 
to calculate the differential 

192 V)* 



dV 



9V dx 

dx 2 dx 2 



dV dA 
dx A 
dV 

-ox—dW' 

ox 



dV , 

dV 1 2 2 d 2 V dVi 
dx 2 dx 2 dt 



dt 



dV 
ax—dW Q 
dx 
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where the last equality holds true, because V is Q-martingale. 

Now we suppose that (t, x) e S. The value of the option is defined by 

V(t,x t ) = pe- qt (l-x t ). 

Hence the differential dV = —pqe~~ qt (l — x)dt — pe~ qt dx has the form 

~ dA 
dV = -pe~ qt x— + pe~ qt {rx - q)dt + pe~ qt axdW Q . 



For both regions we have the following equation 
dV{t,x t ) = -pe- qt l s {t,x t )(x t ^ 



[rx t 



q)dt) + dMj 2 , 



(6) 



where M t is a Q- martingale. Integrating ([6l) from t to T and taking expectation we have 



E; 



V(T,x T ) -V(t,x t ) 



-E; 



1 'dA u 



pe qu x u l s (u,x u ) 



q 



T 



dM? 



V(t,x t ) = E« e- qT (p(l-x T ) 



E; 



pe qu x u l s {u } x u ) (—^ - (r - — )du 

\ -/l 7 , Xn, 



=v(t,x t ) 



=e(t,x t ) 



which completes the proof of Theorem 12.11 



□ 



It is worthwhile noting that the above expression for the value V(t, x t ) = V am (t, x t ) of 
an American-style Asian option can be restated as follows: 



V am {t, x) = V eu (t, x) + E t s / l s (u, x u )f b (u, x u )du 



where V eu (t, x) = E® e qT (p(l — xt 



stands for the price of the European-style Asian 

option and the term l s (u,x u )fb(u,x u ),u G [0, T], represents a surplus bonus for the 
difference between American and European-style of the Asian option. 



3 Early exercise behavior close to expiry 

3.1 Limit of the early exercise boundary at expiry 

In this section we determine the position of the early exercise boundary x* T at expiry 
T. The result is stated for a wide class of integral equations for pricing American-style 
options. 
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Theorem 3.1 Consider an American- style (call or put) option V am on the underlying y 
with the stopping and continuation regions defined by the sets S ^ 7 and C ^ 0, resp. 
Let yl = dS(t, .) = dC(t, .) for t G [0, T] be the early exercise boundary function. Suppose 
that the value of V am is given by the equation 

V am {t,y t ) = V eu {t,y t ) +E t / l s (u,y u )f b (u,y u ) du , (7) 



where V eu denotes the price of the corresponding European- style option and f b (t,y) is a 
continuous function representing the early exercise bonus such that the equation fb{T, y) = 
has a unique root y* > 0. Furthermore we suppose that 

V am (t, y) > V am (T, y) = V eu (T, y) > 0, V am (t, y) > V eu (t, y) for any t G [0, T], y > 0, 

and the function [0, oo) 3 y h- > ^^(T,y) G R is continuous except of the set ATM = 
d{y > 0, V eu (T, y) > 0}. Then only one of the following cases can occur: 

1. y* T G ITM = {y>0, V eu (T } y) > 0}. In this case f b (T,y^) = 0. 

2. y* T G ATM = d{y > 0, V eu (T, y) > 0}. In this case f b {T, y^)>0. 

Remark 3.1 The abbreviation ITM stands for the so-called in-the-money set whereas 
ATM (the boundary of ITM) denotes the so-called at-the-money set. We denote by OTM 
the out-the-money set - the complement of sets ITM and ATM, i.e. OTM = (ITM U 
ATM) . 



Proof: [of Theorem 13.1] We have 



T-t 



ls(u,y u )fb(u,y u ) du = j^—^ (V am (t,y t ) - V eu (t,y t )) > 0, 



for any t G [0,T). In the limit t — > T, we can omit the conditioned expected value 
operator K t and we obtain lg(T,y T )f h (T,y T ) > 0. Since (T, y^,) G dS and the function 
/b is continuous we obtain f h (T,y^) > 0. 

Part 1). Suppose that y* T G ITM. We shall prove f h (T,y^) = 0. Notice that in 
the stopping region S we have the identity V am (t,y) = V am (T,y) for any (t,y) G S and, 
consequently, ^^(t, y) = 0. Take any (T,y T ) G S. Similarly as in above, in the limit 
t — > T we can omit the conditioned expected value operator E t to obtain 

= ^(T,y T ) = ^(T,y T ) ~ fb(T,y T ). 

Since the function y \—> (T, y) is assumed to be continuous in ITM and y^ G ITM, 
we have ^(T j2 £) = f h (T,y^). 

On the other hand, take any (T, y T ) G C. Since V am (t, y) > V am (T, y) for any < t < T 
we obtain 9 g° m (T, y^) < 0. Therefore, in the limit t — > T, we have 

> ^(T,y T ) = ^(T,y T ) - l S (T,y T )f h (T,y T ) = ^(T,y T ) 
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because (T,y T ) G C. Due to continuity of ^p-(T, y) for y G JTM we finally obtain 
> ^Jf (T,?/*) = f h (T,y T ). Hence f h {T,y T ) = and the proof of Part 1) follows. 

Part 2). It suffices to prove that the case y T G OTM cannot occur. Suppose to the 
contrary y T G OTM. Since the set OTM is open we can argue similarly as in the first 
part of the proof to obtain /b(T, Ut) = 0- Since we have assumed uniqueness of the root 
of the equation f^(T,y) = then, in some neighborhood of y T there exists yr G OTM 
such that fb(T,y T ) > and (T,y T ) G S. We have 

dV am / r-w-\ - — . \ dV f 







■(T,y T ) 



(T,Vt) ~ fb(T,y T )- 



dt K dt 
Thus ^{T,y T ) > 0. For any y G OTM we have V eu (T,y) 



0. Hence V eu (T,y T ) = 
and V eu (t, yr) > for all < t < T, a contradiction. Thus y T G" OTM and the proof of 
theorem follows. □ 



As a consequence of Theorem 13.11 we obtain the starting position of the early exer- 
cise boundary for American-style of Asian options with various types of the strike price 
averaging method. 

Corollary 3.1 The value of the limit of early exercise boundary at expiry for the floating 
strike Asian option is summarized in the Tabled In the case of a geometric averaging, it 
follows from (TJP that fb(T, xt) = e~ qT ( — lnxy — rxr + q) such that x T = xt G ITM 
is a solution of the transcendent equation 



. „ qT 
InxT = 



rT. 



The formula for the limit of the early exercise boundary at the expiry ([8]) for geometric 
averaging is the same as presented by ( Wu et al. 1999) and (Detemple 2006, p. 69). Notice 
that the same values of the limit of early exercise boundary at expiry for the continuous 
arithmetic average type of an Asian option is derived also in (Dai and Kwok 2006). 



Table 3: The limit of the early exercise boundary position x T at expiry t = T (xt solves ([5])). 





arithmetic average 


geometric average 


weighted arithmetic average 


put 


max i , 1 


max (iy, 


maX ^(l-e-Wj+A'V 


call 


mm (- , 1 


min (xt, 1J 


min (l(l-^ xt ) + \ A 



3.2 Integral equation for pricing Asian options 

In this section we calculate the approximate formula for the American-style Asian option 
with various floating strike averages. The next lemma will be useful in calculations to 
follow. 

Lemma 3.1 (Hansen and J0rgensen 2000) Let to = \nQ ~ Af(a,j3 2 ) and define 7 = 
a+^-injf ^ where K>Q We have 

E [ l {pn>pK } ] = *(p(7 -/?)), E[l {pn > pK} mnn] = e a+ ^((« + /5 2 )$(p 7 ) +pP$(l)), 
E [l{pn> P i<:}^] = e a+ ^$(p7), - K)) + ] = p(^e a+ ^~^(p^) - K$>(p{j - 
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Figure 1: Isolines of the the limit of the early exercise boundary at expiry T of call option 
for the continuous arithmetic (left) and the continuous geometric average (right). 



where p G { — 1, 1} and $( 
density functions, resp. 



and $(■) are standard normal cumulative distribution and 



3.2.1 Geometric average 

In this section we recall the integral equation for pricing American-style of Asian geomet- 
rically averaged floating strike options. It was derived for the case q = by (Hansen and 
J0rgensen 2000) and for the general case q > by (Wu et al. 1999). 

A? 

Lemma 3.2 (Wu et al. 1999) In the case of geometric averaging, the variable x\ = has 
log-normal (conditioned) distribution In x^Tt ~ Af(a t , u , /3t, u ) > where u > t and parameters 
ott,u = a(t,u,x t ) and j3 t ,u = (3(t,u) are defined by 



a 9 t,u = -lnxf 
u 



u 



t 2 



2u 



ff+y), Pin 



a 



V3 



t 3 . 



u 



Now one can apply Lemma 13.11 in order to calculate the formula for option with the 
geometric averaging. Recall that for the floating strike Asian call or put option, the 
stopping region S = {{t,x),x > 0,px^ > px}, where x* t is the exercise boundary and 
p = 1 for the case of a call option whereas p = — 1 for a put option. If we insert the 
expression -jf for the geometric average (see Table [1]) into (0J and ([5]) we obtain the 
formula for the European-style option 



v 9 (t,x t ) 



E 



pe 



-qT 



(p(i - A\ 



-qT, 



a 



t.r 



a 9 - 
e a t,T 



a 



t.r 



+ A 9 
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and the value of the American early exercise bonus premium 

cT 



r f /dA 9 \ 

e 9 (t,x t ) = E?[J pe-« u xll s {u,xi)(^ - {r - q{xl)- l )du) 

= jf pe-*>E* [l {px * u>pxi} ( - \xl In xi - rx 9 u + q 



T 



pe-^[q*Wlu-l!,«)) 



f - ( 



where $(■) is the standard normal cumulative distribution function and 

g a 9 t tU -lnxl 

1t,u ng ' Pt,u m 

Pt,u 

Returning to the original variables we obtain the formula of American- style floating 
strike Asian option with geometrically averaged floating strike: 

V 9 {t, S, A) = Se qt V 9 {t, A/S) = Se qt {v 9 {t, A/S) + e 9 {t, A/S)). 

If we formally set value of the continuous dividend rate to zero, i.e. q = 0, the result is 
identical to the expression obtained in the paper (Hansen and J0rgensen 2000)). 

3.2.2 Approximation for the arithmetic average 

Unfortunately, in the case of an arithmetically averaged floating strike Asian option the 
probabilistic distribution function of the arithmetic average cannot be expressed in an 
explicit way. Following (Hansen and J0rgensen 2000) we approximate the probabilistic 
distribution of the variable = for the continuous arithmetic average A® by the 
log-normal conditioned distribution, i.e. lna;*!^ ~ J\f(a% u , {Ptu) 2 ) a ^ time t, where 



< u = 21nEf [<] -~lnE t s [«) 2 ], /3« u = ^/ln E t s [(a-)*] - 2 In E t s . (9) 
Lemma 3.3 Consider the variable x u = 4^, where A u and S u are defined as the arith- 

u 

metic average (see Table{l\) and as in resp. First two conditioned moments E®[x^] 
and Ep[(x°) 2 ] of x u entering the expressions for the functions ct1 u = a a (t,u, xf) and 
Ptu = ^ a {^-> u -> x t) can be calculated, fort < u, as follows: 

E, 2 r<l = tfte-tr-M*-*') + - X — (l - e ~(r-g){u-t)\ (10) 
L J u [r — q)u 

T- 2 2 0+f,-(r-q)(u-t) 
E t 2 r«) 2 l = ( X ?) 2 L- e ~2(r-<l-^)(u-t) +x a^_ / x _ e -(r-,)(„-m (U) 

u u 2 (r — q) 

(r - q - a 2 ) - 2(r - q - ^-) e -^-^ u ~^ + (r - q^-^-i-^-t) 
u 2 (r — q)(r — q — ^-){r — q — a 2 ) 
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Proof: Following the lines of the derivation of E^[a;^] from (Hansen and J0rgensen 2000) 
adopted for a general dividend yield q > we obtain the first moment 



E 



i2 



x: 



x+ — e 
u 



-(r-q)(u-t) 



_ e -(r-q)(u-t) 

(r — q)u 



Although we follow the proof by (Hansen and J0rgensen 2000), we have to make a slight 
correction in the derivation of the second conditioned moment. Using the definition of St, 
we have, for all v G [t,u], 



E? 


S v 






-S u - 





Ar-q+^Kv-u)+a(WZ-W£) 



,{r-q)(v-u) 



We need to simplify the expression for the the second conditioned moment 



E 



Q 



x 



a\2 



U z 



= E? 

St St 
S u S u - 



r/i rs 



u Jo S u 
t 
u 



-dv 



2-i 



■ 2x t — E t s 



St 

Su- 



E; 



S v 
-S u - 



dv + 



u 



E 



S z S v 
S u S u - 



dv dz. 



Assuming that u > z,v and let m = min {z, v} and M = max {z, v}, we have 

S z s. 



«T(2(wa-wjg)+(wjg-wO)) 



E; 



S z S v 
S u S u 



j(r-q)(z+v-2u)+a 2 (u-M) 



We have calculated all expressions we need to evaluate the second conditioned moment. 
If we put all together and perform necessary calculation we obtain (lllj) and the proof of 
lemma follows. □ 



Remark 3.2 If we formally set the value of the continuous dividend rate q = in 
Lemma l3.3\ we obtain almost identical expression to that of (Hansen and J0rgensen 2000) 

except of the second moment Ep (x") 2 entering (TJ|). The expression 

m«) 2 ] H j = (s?) 2 4e- 2 (-g)(«-'> + xt 2t 2 f (U ~l (1 - e-(- 2 >^>) 

u u z (r — cr z ) 

T 2 



(T 2 ) - 2(r - 4)e~ r( "-* ) + re-^-^X"-*) 



u 2 r(r — ^-)(r — a 2 ) 

6y (Hansen and J0rgensen 2000) differs from our (Tl\) in the second summand where the 
term r — a 2 is replaced by r in both in the denominator and exponent. 

Now we can return to the problem of valuation of an Asian option. First we replace 
the general form of the average in (J5J) by the expression for the arithmetic average defined 
as arithmetic average in Table [TJ The stopping region S is the same as for the case of 
geometric averaging. 
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Using Lemma 13.11 we calculate the value of both and (jSJ) • The European part of 
the option has value 



v a (t,x) 



E? 



- 9T 'p(l -x a T ) 



pe 



and the American early exercise bonus premium 

r f T /dA a 
?(t,x) = Ef / p e-V u l s (u,x a a )(-^-(r-q(x a a )- l )du 



^ Jt 



Aa 



pe-^Ef l {pxl>pxi} -(1 - <) - rxl + q 



u 



du 



pe 



- qu 



where $(•) is the standard normal cumulative distribution function and 

a?„ - In 

o a ^ Pt,u- 

Pt,u 

Returning to the original variables we have the approximate value of American-style 
Asian option with a continuous arithmetic averaging 



V a {t, S, A) = Se qt V a {t, x) = Se qt {v a {t, x) + e a {t, x)). 



3.3 Expansion of the early exercise boundary close to expiry 

Throughout this section we shall assume the structural assumption on the interest and 
dividend rates: 

r>g>0. (12) 

We shall calculate an approximation of the call option early exercise boundary function for 
a call option. The approximation is obtained by the first order Taylor series expansion in 
the a/t variable, where r = T — t is the time to expiry. To approximate the early exercise 
boundary function by the Taylor expansion, we need to calculate the first derivative of x^ 
at expiry T with respect to \/T — t variable. Following (Kuske and Keller 1998, Dewynne, 
Howison, Rupf and Wilmott 1993, Sevcovic 2001) we propose an approximation of the 
early exercise boundary x* t in the form 

x* = x* T (l + Cy/T - t) + 0{T -t) as t ^ T, 

where C G K. is a constant. To calculate C, we use the condition of smoothness of the 
early exercise boundary of the call option across the early exercise boundary - smooth 
pasting principle (c.f. (Kwok 1998, Dai and Kwok 2006)). Since V(T, x) = e~ qT (l — x) + 
we have 



13 



BV Bv rfe f T 

-l = e*—(t,xt) = e«^(t,xt)+e«-^(t,x* t )=v x (t,x* t )+ J e Ix (t, x* t , u, <) du. (13) 

In further derivation we use following limits (according to the Lemma 13.21 Lemma 13.31 
and Table E]). We recall that a t)U = a(t,u,x), = f3 a (t,u,x a ) and (3 9 t)U = (3 9 (t,u). 

lim cut-t t = h m oct-t t-t(i-9) = «tt = hix^ < 0, 
lim Pt-t,t = lim (3 T -t,t-t{i-6) = Pt,t = + , 



OT-rT 



lim >, /3r '" T/ = 0, (14) 
lim lnx ^^ + /l(T V /r ( 1 ~ 9 )) - ar-r,T-r(i-e) , 1 - Vj - g 

#T-t,T— r(l-0) \/~0~ 



lim d x a T _ T)T = lim d x a T - T ,t-t(i-o) = ~ 



lim (3t-t.t d T (3 T -T,T = tt ■ 

Since we have assumed r > g > we have < iy < 1 (see Table [3]). Notice that 
both a and /3 have polynomial order in r and the derivative of the normal cumulative 
distribution function (i.e. the probability density function) has exponential order in r 
variable. In both derivations we have used several properties of the derivative of normal 
cumulative distribution function $(#), e.g. $'(x) = <&'(— x), <&"(x) = —x&(x) and 

*'(§ + &) = e- a -^$'(f). 

The following lemma will be useful in derivation of asymptotic behavior of the early 
exercise. Its proof is straightforward and follows from monotonicity of the right-hand side 
of equation fllSH in the h variable. 



Lemma 3.4 The implicit equation 

= 1 - f •( " h'-^dO + h I' - h'-^dO (15) 



j o \/£ ' Jo Ve ^ Ve > 

has the unique solution h* having its approximate value h* = —0.638833. 

Notice that the first order asymptotic expansion as t — > T of the early exercise 
boundary Sf(t) ~ S/(T)(1 + 0.638833 o~\/T — t) for the American call option derived 
by (Dewynne et al. 1993, Sevcovic 2001) contains the same constant —h* = 0.638833 
where h* is a solution of (I15D. 



Remark 3.3 For the early exercise boundary function x* = xl(T,r,q,a 2 ) as a function 
of the model parameters T,r,o~ 2 > 0, q > 0, we have the following scaling property: 



x* t (T,r,q,a 2 ) = x%{1, rT, qT, a 2 T). 
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3.3.1 Geometric average 

We calculate the derivative of the European part of the expression (for the call option). 
Recall that a depends on the variable x, but f3 does not depend on this variable. 

Now we calculate the derivative of the integral function of American-style option bonus: 

~* u *\ -,(«-*) 1 ( {-qv. + x* u ru + x* u kix* u + x*M >u ) 2 \ Anx* u -a 9 t y 
& Ix (t,x,u,x u ) = e q( >-U —g 1*1 

^(™+< + (Ag 2 + i)*( 



— e 




lim d T a 9 T _ rT d x a 9 T _ rT + d\ T a 9 T _ rT = i (q 



We want to determine the behavior of the early exercise boundary near the expiry 
T. The limit in the expression ( fl3l) leads to the trivial identity. By rearranging all of 
elements on the right side of the equation, we have an expression of order T — t. We 
substitute r = T — t and x* T = G and divide the equation by r. We have: 

i + e, ( r-T,G(i + w?)) = lim _ G(1 + , !ffV7)) . 

According to the Lemma [3.21 and Table [3] we have: 

1 _ I _ ^ 

GV G~T~Y 

The only non-zero elements of the first partial limit are the elements multiplied by the 
cumulative distribution function. 

r->o T T->odr\ dx J G T 

The second term represents the limit of the integral part divided by r. If we substitute 
u = T — t(1 — 0), then we obtain 

/t-t - ^ ^l 1 + ^v 7 ^)), u, G{1 + haVT=u))du 
lim 

= / limef^T - r,G(l + hay/r)),T - t(\ - 6),G{1 + ha^?VT^d))dd. 
Jo 



The last expression then can be simplified, using limits ffl4|) . equation (jo]) for the limit of 
early exercise boundary at expiry G, i.e. r = ^ — ^P, and by calculating the limit (using 
L'Hospital rule) of the expression multiplied by the derivative of the cumulative density 
function. The final limit of the integral has form 



Hmef^T - r, G(l + ha^)),T - r(l - 9), G(l + ha^Vl^O)) 




Integrating ([15]) with respect to 6 G [0, 1], putting both partial limits together, dividing 
by the nonzero constant q + if and by Lemma 13.41 and Table [3], we finally obtain 

x * = G(l + h*aVT - t) + 0(T -t) as t -> T, where h* = -0.638833. 
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3.3.2 Arithmetic average 

The derivation for the case of the arithmetic average is very similar to the geometric one. 
We calculate the derivatives of both parts of the value function. The European part 

d~^^ -a(T-«^( <T\dPtT 



6 ftp ^ T )\ dx + ^ dx 

and the American-style bonus 




,fl nx* u -al u \ dffi tU 

PL / ^ 



1\ „ , tft%) 2 /lnx*-a?„ \ /daf„ d8?„* 



The rest of derivation was performed following the same steps used in section 13.3. 1[ 
The main difference in the derivation is that for the arithmetic average also the expression 
Ptu = P a (t, u, x) depends on the variable x. Thus, according to the Lemma [3731 and Table[3] 
we need to calculate following limits 



lim d T a a T _ T ^ T d x a^_ TtT + d\ T a a T _ rT = —(q- r 



I / 1 _ o2 

f ~ Y 



limd x (%_ T>T = hmc^_ TiT _ T(1 _ e) = limd x f%_ TtT d x ft_ TiT + /^_ t , t <%^t-.,t = 0. 

We recall that for the case of a continuous arithmetic average we have x*(T) — A = 

Since we have assumed r > q > we obtain A < 1 (see Table [3]). 

The derivation leads to the same equation as in Lemma 13.41 when multiplied by the 
constant r + ^. In summary, we obtain the following approximation of the limiting 
behavior of the early exercise boundary near expiry also for the arithmetic average. 

x*(t) = 1 + qT (l + h*aVf -t) + OCT - t), ast^T, where h* = -0.638833. (17) 
1 + rT 



4 Transformation method for Asian call options 

The purpose of this section is to propose an efficient numerical algorithm for determin- 
ing the free boundary position x\ for American-style of Asian options. Construction of 
the algorithm is based on a solution to a nonlocal parabolic partial differential equation 
(PDE). The governing PDE is constructed for a transformed variable representing the 
so-called 5-synthetised portfolio. Furthermore, we employ a front fixing method (refereed 
also to as Landau's fixed domain transformation) developed by (Stamicar, Sevcovic and 
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Chadam 1999, Sevcovic 2001) for plain vanilla options as well as for a class of nonlin- 
ear Black-Scholes equations (Sevcovic 2007, Sevcovic 2009). At the end of this section 
we present numerical results and comparisons achieved by these methods to the recent 
method developed by (Dai and Kwok 2006). 

First, we recall the partial differential equation for pricing Asian options (cf. (Kwok 
1998)). We assume the asset price dynamics follows a geometric Brownian with a drift g, 
continuous dividend yield q > and volatility a, i.e. dS = (g — q)Sdt + aSdW where W 
is the standard Wiener process. If we apply Ito's formula to the function V = V(t, S, A) 
we obtain 

n . ,'dV cr 2 n2 d 2 V\ , dV , n dV, A 

dV=[— + —S 2 —-)dt + —dS + —dA. 18 



dt 2 dS 2 J dS OA 

Recall that for arithmetic, geometric or weighted arithmetic averaging we have dA/A 
f(A/S,t)dt where the function / = f(x,t) is defined as follows (see Tab. [2]): 

arithmetic averaging, 



f(x,t) = { Hf 2 geometric averaging, (19) 

exponentially weighted arithmetic averaging. 



A(x~ 1 — 1) 
l_ e -xt 



Inserting the expression dA = Af(A/S,t)dt into (ITS]) and following standard arguments 
from the Black-Scholes theory we obtain the governing equation for pricing Asian option 
with averaging given by (IT9l in the form: 

dV a 2 n d 2 V dV dV 

where < t < T, S,A> (see e.g. (Dai and Kwok 2006, Kwok 1998)). For the Asian 
call option the above equation is subject to the terminal pay-off condition V(T, S, A) = 
max(S" — A,0), S,A > 0. It is well known (see e.g. (Kwok 1998, Dai and Kwok 
2006)) that for Asian options with floating strike we can perform dimension reduction by 
introducing the following similarity variable: 

x = A/S, W(x,r) = V(t,S,A)/A 

where r = T — t. It is straightforward to verify that V(t, S,A) = A W(A/S, T — t) is a 
solution of (1201 iff W = W(x, r) is a solution to the following parabolic PDE: 

dW a 2 d ( 2 dW\ , . dW p . m f dW\ . . 

x 2 —- + (r - q)x— f(x, T-t)[W + x—- + rW = 0, (21) 



dr 2 dx \ dx J dx \ dx 

where x > and < r < T. The initial condition for W immediately follows from the 
terminal pay-off diagram for the call option, i.e. W(x, 0) = max(x~ 1 — 1,0). 



4.1 American-style of Asian call options 

Following (Dai and Kwok 2006) the set £ = {(*, S, A) e [0, T] x [0, oo) x [0, oo), V(t, S, A) = 
V(T, S, A)} is the exercise region for American-style of Asian call options. In the case 
of a call option this region can be described by the early exercise boundary function 
S f = S f (t,A) such that £ = {(t,S,A) E [0,T] x [0,oo) x [0,oo), S > S f (t,A)}. For 
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American-style of an Asian call option we have to impose a homogeneous Dirichlet bound- 
ary condition V(t, 0, A) = 0. According to (Dai and Kwok 2006) the C 1 continuity con- 
dition at the point (t, Sf(t, A), A) of a contact of a solution V with its pay-off diagram 
implies the following boundary condition at the free boundary position Sf(t, A): 

?^(t,S f (t,A),A) = l, V(t,Sf(t,A),A) = S f (t,A)-A, (22) 

for any A > and < t < T. It is important to emphasize that the free boundary 
function 5/ can be also reduced to a function of one variable by introducing a new state 
function ctS follows: 

S f (t,A) = A/x* t . 

The function t ^ x+ is ci free boundary function for the transformed state variable x = 
A/S. For American-style of Asian call options the spatial domain for the reduced equation 
( 121]) is given by 1/ q{t) < x < oo, r G (0, T), where g(r) = l/x^_ T . Taking into account 
boundary conditions (T22]) for the option price V we end up with corresponding boundary 
conditions for the function W: 

1 dW 1 1 

W(oo,t) = 0, W(x,t) = - - 1, -z-(x,t) = ; at x = — - (23) 

x ox x 2 q{t) 

for any < r < T and the initial condition 

W(x, 0) = max(a; _1 -1,0) for any x > 0. (24) 



4.2 Fixed domain transformation 

In order to apply the Landau fixed domain transformation for the free boundary problem 
( f2~Tl) . ( |23i) . ( 1241) we introduce a new state variable £ and an auxiliary function II = n(£, r) 
representing a synthetic portfolio. They are defined as follows: 

dW 

£ = In (f?(r)x) , n(£, r) = W(x, r) + ar-^fo r) . 

Clearly, x G (^(r) _1 , cxd) iff £ G (0, cxd) for r G (0, T). The value £ = oo of the transformed 
variable corresponds to the value x = oo, i.e. 5 = when expressed in the original 
variable. On the other hand, the value £ = corresponds to the free boundary position 
x = x%, i.e. 5 = Sf(t, A). After straightforward calculations we conclude that the function 
II = II(£, r) is a solution to the following parabolic PDE: 

where the term a(£, r) depends on the free boundary position g. The terms a, b are given 
given by 



«(M = 44+r-g-^-/(e«Mr),T-r), 6(£, r) = r-^- (x/(x, T - r)) 
g(r) 2 ax 



. (25) 

- e ^ 
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Notice that &(£, r) = r + 1/(T — r) in the case of arithmetic averaging, i.e. f(x,t) 
The initial condition for the solution II can be determined from (EH) 



n(£,o) 



-1 £< ln2(0), 
£>lng(0). 



Since d x W(x, r) = — ^ and W(x, r) = | — 1 for x = and VT(oo, r) = we conclude 
the Dirichlet boundary conditions for the transformed function Il(£, r) 

n(0,r) = -l, n(oo,r) = 0. 

It remains to determine an algebraic constraint between the free boundary function g(r) 
and the solution II. Similarly as in the case of a linear or nonlinear Black-Scholes equation 
(cf. (Sevcovic 2007)) we obtain, by differentiation the condition W(-^,t) = g(r) — 1 
with respect to r, the following identity: 

dg dW , , , i . , , ,_ 2 n dg . . dW , . . , . 
(r) = -^(g(r) 1 ,r)(- f? (r) 2 ) -£(r) + -jt-(q(t) 1 ,r). 



dr dx ' dr d 



T 



Since c^W^^r) 1 ,r) = — f?(r) 2 we have Wr{x,r) = at x = g(r). Assuming continuity 
of the function IT(£, r) and its derivative II^(^, r) up to the boundary £ = we obtain 

2 d 2 W f dU /n , , . . . . , . , 

x -^{ X ,V ^> —{0,t) +2g{T), x—(x,t) -»• -g(r) as ac -> g(r) . 

Passing to the limit x — > f?(r) _1 in (12T1) we end up with the algebraic equation 

qg(r) - r + f(g(r)~\ T - r) = £25(0, r) (26) 

for the free boundary position ^(r) where r e (0, T]. Notice that, in the case of arithmetic 
averaging where /(p(r) _1 , T — r) = (g(r) — 1)/(T — r), we can derive an explicit expression 
for the free boundary position g(r) 

l + r(T-r) + 4(T-r)f (0,r) 
= l + g (T-r) ' 0<T<T ' 

as a function of the derivative ^11(0, r) evaluated at £ = 0. The value g(0) can be deduced 
from Theorem 13. II For the arithmetic averaging we have (see also (Dai and Kwok 2006)) 
the following expression: 

(1+rT 

g(0) = max 

In summary, we derived the following nonlocal parabolic equation for the synthesized 
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portfolio n(£, r) 



<9IT /(h ,011 



a 2 d 2 U 



2 d£ 2 

with an algebraic constraint 



+ 6(£,t)II = 0, 0<r <T, £>0, 



qg(r)-r + f{ e (T)- l ,T-T) 



y^ (0 ' r) 



n(£,o) 



< r < T, 



subject to the boundary and initial conditions 

n(0,r) = -l, n(oo,r)=0, 

-1 for £ < ln(e(0)), 
for£>lnte(0)), 

where a(£, r) and &(£, r) are given by (|25|) . 

and the starting point g(0) — is given by Theorem 13.11 



(27) 



4.2.1 An equivalent form of the equation for the free boundary 

Although equation ( l26i) provides an algebraic formula for the free boundary position q{t) 
in terms of the derivative ^n(0, r) such an expression is not quite suitable for construction 
of a robust numerical approximation scheme. The reason is that any small inaccuracy in 
approximation of the value ^11(0, r) is transferred in to the entire computational domain 
£ G (0, oo) making thus a numerical scheme very sensitive to the value of the derivative 
of a solution evaluated in one point £ = 0. In what follows, we present an equivalent 
equation for the free boundary position q{t) which is more robust from the numerical 
approximation point of view. 

Integrating the governing equation (j2Tj) with respect to £ G (0, oo) yields 



d 



IM£ 



o~ 



d 2 U 



dc+ / 6(£,r)IM£ = 0. 



drJo Jo " 2J «9£ 2 

Now, taking into account the boundary conditions 11(0, r) = — l,II(oo,r) = 0, and con- 
sequently c^n(oo, t) = we obtain, by applying condition (|2"6T) . the following differential 
equation: 



d ( f°° \ <j f 

— Mn q(t) + J II(£, r)dn + 9e (r) - q - — + 

In the case of arithmetic averaging where f(x,t) = 



r- f 



-i 



d 

dr 



In q(t) 



a 



U(^T)dn+qg(r)-q- — + 



l)/t we obtain 

e(r)e _€ - 1 
T-r 



n(e,r)de=o. 



n(£,r)d£ = 0. 

(2s; 



4.3 A numerical approximation operator splitting scheme 

Our numerical approximation scheme is based on a solution to the transformed system 
(f27|) . For the sake of simplicity, the scheme will be derived for the case of arithmetically 
averaged Asian call option. Derivation of the scheme for geometric or weighted arithmetic 
averaging is similar and therefore omitted. 
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We restrict the spatial domain £ E (0, oo) to a finite interval of values £ E (0, L) where 
L > is sufficiently large. For practical purposes it sufficient to take L pa 2. Let > 
denote by the time step, k = T/m, and, by h = L/n > the spatial step. Here m, n E N 
denote the number of time and space discretization steps, resp. We denote by IP = IP (£) 
the time discretization of n(£,Tj) and gi ~ Q{ T j) where Tj = jk. By H\ we shall denote 
the full space-time approximation for the value II(£j, Tj). Then for the Euler backward in 
time finite difference approximation of equation (|27|) we have 

IP -IP'" 1 ,dW fa 2 p>e-t-l\dW a 2 d 2 W ( 1 . 

+ J—r - — + ^ ~ -r^T- + [r + - IF = 



k <9£ \2 T-Tj J d£ 2 d 2 i \ T - r. 

where c 7 is an approximation of the value c[jj) where the c(r) = |£4 + r — q. The 

solution W = IP'(x) is subject to Dirichlet boundary conditions at £ = and £ = L. We 
set n°(£) = I1(£,0) (see (|2"T)) ). In what follows, we make use of the time step operator 
splitting method. We split the above problem into a convection part and a diffusive part 
by introducing an auxiliary intermediate step IP~2: 
( Convective part) 

Ip-| - IP'- 1 • i 

- + d>d x W "5 = 0, (29) 



(Diffusive pari) 



W -Jp-k (a 2 e>e-t-l\dW a 2 d 2 W , , , 

-7T + ^ -zr - -tt^tt +[r + IF = 0. (30) 



k \ 2 T — Tj J <9£ 2 <9 2 £ V T - r i 

Similarly as in (Sevcovic 2007) we shall approximate the convective part by the explicit 
solution to the transport equation d T Yl + c(r)d(U. = for £ > and r E (tj_i,Tj] 
subject to the boundary condition 11(0, r) = —1 and the initial condition n(£, 7^-_i) = 
rF _1 (£). It is known that the free boundary function q(t) need not be monotonically 
increasing (see e.g. (Dai and Kwok 2006, Sevcovic 2009) or (Hansen and J0rgensen 
2000)). Therefore depending whether the value of c(r) is positive or negative the boundary 
condition n(0, r) = — 1 at £ = is either in-flowing (c(r) > 0) or out-flowing (c(r) < 0). 
Hence the boundary condition n(0,r) = —1 can be prescribed only if c(jj) > 0. Let us 
denote by C(r) the primitive function to c(r), i.e. C(r) = In q(t) + (r — q)r. Solving the 
transport equation d T H + c(r)d^U. = for r E [rj-.i,Tj] subject to the initial condition 
m,r^j) = IP" 1 ® we obtain: n(£,r) = W^-C{r)+C{Tj^)) ii ^ - C (r) + C (r 3 ^) > 

and n(£, r) = —1 otherwise. Hence the full time-space approximation of the half-step 

3 -_i 

solution 2 can be obtained from the formula 

n J "5 = J 11 ^ 1 ^)' if Vi = & ' M + ln^" 1 - (r - q)k > 0, 
* |^ — 1, otherwise. 

In order to compute the value IP' -1 (?7i) we make use of a linear interpolation between 
discrete values n^ _1 , i = 0, 1, n. 

Using central finite differences for approximation of the derivative d x W we can ap- 
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proximate the diffusive part of a solution of (130]) as follows: 

n? - ii- • / i 



k V T-Tj 



n 



.i 



± , g*e-<'-l \ g+i - g-i _ a 2 nj +1 - 211? + nj_x = 
2 T-Tj J 2h 2 //- 



Therefore the vector of discrete values IP = {1^1, i = l,2,...,n} at the time level j G 
{1,2,..., m} is a solution of a tridiagonal system of linear equations 



c^UU + fini + 7 M +1 = \ for i = 1, 2, n, where (32) 

= 1 + (r + fc - (aj + 7f) • (33) 

The initial and boundary conditions at r = and x = 0, L, can be approximated as 
follows: 



n° 



-1 for & < \n((l + rT)/(l + qT)), 
for & >ln((l + rT)/(l + gT)), 

for i = 0,1, n, and IIq = —1, lP n = for j = 1, m. 

Finally, we employ the differential equation (|28|) to determine the free boundary posi- 
tion g. Taking the Euler finite difference approximation of (in g + J °° IT<i£) we obtain 

(Algebraic part) 



In ^ = In 1 + ^(n^ 1 ) - I (W) +k(q+^- qg^ 1 - h(g j -\ IP)) 



(34) 



where lo(II) stands for numerical trapezoid quadrature of the integral J °° II(£)d£ whereas 
Ji(^ _1 , II) is a trapezoid quadrature of the integral / °° — ^-j^— - j II(£)cf£. 
We formally rewrite discrete equations (I3TI) . (|32|) and (I34p in the operator form: 

g> = .F(IF), IP~i = Tig'), A{g j )W = IF"*, (35) 

where lnjF(IP) is the right-hand side of equation f[3~4"|) . T(g>>) is the transport equation 
solver given by the right-hand side of (13T!) and ^4. = A(g>) is a tridiagonal matrix with 
coefficients given by (!33l) . The system (1351) can be approximately solved by means of 
successive iterations procedure. We define, for j > 1, IP' = IF , g^'° = g 3 . Then the 
(p + l)-th approximation of IP and g> is obtained as a solution to the system: 

IP-!,p+i = T{g j ' p+1 ), (36) 

^(^'• p+1 )n J ' p+1 = iF-5>p +1 . 

Supposing the sequence of approximate discretized solutions {(IP' P , g^' p )}^Li converges to 
the limiting value (IP' 00 , ^ 7,0 °) as p — > oo then this limit is a solution to a nonlinear system 
of equations (135]) at the time level j and we can proceed by computing the approximate 
solution in the next time level j + 1. 
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Figure 2: The function q(t) (left). A comparison of the free boundary position x[ = l/g(T — t) 
(right) obtained by our method (solid curve) and that of the projected successive over relaxation 
algorithm by Dai and Kwok (dashed curve). 




Figure 3: A 3D plot (left) and contour plot (right) of the function II(£, r). Profiles of the 
function i1(£,t) for various times r G [0, T]. 

4.4 Computational examples of the free boundary approxima- 
tion 

Finally we present several computational examples of application of the numerical approx- 
imation scheme (l36il for the solution II(£, r) and the free boundary position q(t) of ( l27j) . 
We consider American-style of Asian arithmetically averaged floating strike call options. 

In Fig. [2] we show the behavior of the early exercise boundary function q(t) and the 
function x* t = 1/ g(T—t). In this numerical experiment we chose r = 0.06, q = 0.04, a = 0.2 
and very long expiration time T = 50 years. These parameters correspond to the example 
presented by (Dai and Kwok 2006). As far as other numerical parameters are concerned, 
we chose the mesh of n = 200 spatial grid points and we have chosen the number of time 
steps m = 10 5 in order to achieve very fine time stepping corresponding to 260 minutes 
between consecutive time steps when expressed in the original time scale of the problem. 

In Fig. [3] we can see the behavior of the transformed function II in both 3D as well as 
contour plot perspectives. We also plot the initial condition II(£, 0) and five time steps of 
the function £ i— > Il(£, Tj) for Tj = 0.1, 1, 5, 25, 50. 

A comparison of early exercise boundary profiles with respect to varying interest rates 
r and dividend yields q is shown in Fig. HI A comparison of the free boundary position 
x* t = l/g(T — t) obtained by our method (solid curve) and that of the projected successive 
over relaxation algorithm by (Dai and Kwok 2006) (dashed curve) for different values of 
the interest rate r is shown in Fig. H] (right). The algorithm of (Dai and Kwok 2006) is 
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Figure 4: A comparison of the free boundary position x* for various dividend yield rates 
q = 0.04, 0.03, 0.025 and fixed interest rate r = 0.06 (left). Comparison of x* for various interest 
rates r = 0.06, 0.04, 0.02 and fixed dividend yield q = 0.04. Dots represents the solution obtained 
by Dai and Kwok (right). 



based on a numerical solution to the variational inequality for the function W = W(x, r) 
which a solution to (l2Tj) in the continuation region and it is smoothly pasted to its pay-off 
diagram (|23|) . It is clear that our method and that of (Dai and Kwok 2006) give almost 
the same results. A quantitative comparison of both methods is given in TableHJfor model 
parameters T = 50, o = 0.2, q = 0.04 and various interest rates r = 0.02, 0.04, 0.06. We 
evaluated discrete L°°(0,T) and L 1 (0,T) norms of the difference x * t ,trans — x * t ,psor between 
the numerical solution x* t ,trans ,t e [0,T], obtained by our method and that of (Dai and 
Kwok 2006) denoted by x* t ,psor . We also show the minimal value min tg [ 0T ] x * t ,trans of the 
early exercise boundary. 



Table 4: Comparison of PSOR and our transformation method for T = 50, a = 0.2, q = 0.04. 





r = 0.06 


r = 0.04 


r = 0.02 


|| *. trans *,psori 

iPt ~ x t iioc 


0.09769 


0.03535 


0.05359 


|| *, trans *,psorn 
\\ x t ~ X t 111 


0.00503 


0.00745 


0.01437 


*, trans 

mm x t 


0.52150 


0.57780 


0.63619 



Finally, in Fig. Owe present numerical experiments for shorter expiration times T = 0.7 
and T = 1 (one year) with zero dividend rate q = and r = 0.06, a = 0.2. We also present 
a comparison of the free boundary position x\ = l/g(T—t) and the analytic approximation 
( TT7|) for parameters: r = 0.06, q = 0, a = 0.2 and T = 1. It is clear that the analytic 
approximation (I17p is capable of capturing the behavior of x* only for times t close to 
the expiry T. Moreover, the analytic approximation is a monotone function whereas the 
true early exercise boundary x* t is a decreasing function for small values of t and then it 
becomes increasing (see e.g. Fig. [5]). 



Conclusions 

In this paper we analyzed American-style Asian options with averaged floating strike. 
We focused on arithmetic, geometric and weighted arithmetic averaging of the floating 
strike price. In the first part of the paper we derived an integral representation of the 
call and put option prices and we provided an integral equation for the free boundary 
position. We analyzed the behavior of the early exercise boundary close to expiry. We 
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Figure 5: The free boundary position for expiration times T = 0.7 (left) and T = 1 (center). A 
comparison of the free boundary position with its the analytic approximation (dashed line). 

proposed a general methodology how to determine the early exercise position at expiry. 
We furthermore derived the asymptotic formula for the early exercise boundary close to 
the expiry. The second part of the paper was devoted to construction of a robust numerical 
scheme for finding an approximation of the early exercise boundary. Applying the front 
fixing method, we derived a nonlocal parabolic partial differential for the synthetised 
portfolio and the free boundary position. Using an idea of the operator splitting technique 
we moreover constructed a numerical scheme for numerical solution of the problem. The 
capability of the method has been documented by several computational examples. 
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